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Abstract. The elastic theory of quasicrystals considers, in addition to the "normal" displacement field, 
three "phason" degrees of freedom. We present an approximative solution for the elastic Green's function 
of icosahedral quasicrystals, assuming that the coupling between the phonons and phasons is small. 

PACS. 61. 44. Br Quasicrystals - 62.20.Dc Elasticity, elastic constants 



1 Introduction 

The displacement field of an infinite linear elastic medium 
caused by external forces can be calculated immediately, 
when the elastic Green's function is available. Unfortu- 
nately, analytical expressions for the Green's function are 
known only for special cases. Lord Kelvin [|IJ has pro- 
vided a solution for the isotropic medium, Lifshitz and 
ROSENTSVEIG || and Kroner [|[ for hexagonal crystals. 
Approximations have been presented by Dederichs and 

\ Leibfried f| for cubic crystals and by Kroner || for 
arbitrary anisotropic systems. 

In the case of quasicrystals || , three additional phason 
degrees of freedom lead to a more general theory of elas- 

■ ticity pj-|lO|] . The corresponding Green's functions for pla- 
nar pentagonal (De and Pelcovits decagonal and 
dodecagonal (Ding et al jn]]) quasicrystals are known. 
However, the solution for icosahedral quasicrystals with 
five elastic constants (two for phonons, two for phasons, 
and one for the coupling between phonons and phasons) 

[ has not been established yet. 

It is the intention of this paper to derive an analytical 
approximation of the elastic Green's function for icosahe- 
dral quasicrystals. 

The paper is organised as follows. In section || we recall 
the elastic theory of icosahedral quasicrystals. We present 
the Green's function in the case of the "spherical approxi- 
mation" and for vanishing coupling between phonons and 
phasons. Starting from this solution and using perturba- 
tion theory as well as an invariant basis for the Green's 
function tensor, we derive in section |^ an approximative 
solution for finite phonon-phason interactions. The last 
section presents the phonon and phason fields for some 
special cases. 



2 Elastic theory for icosahedral quasicrystals 
2.1 Equations 

The elastic theory of quasicrystals [H-§ [l2[ [l3| studies, 
apart from the classical "phonon" displacement field u(x), 
also the "phason" displacement field w(x). Icosahedral 
quasicrystals can be described as a three-dimensional cut 
through a six-dimensional cubic crystal. The six-dimen- 
sional space divides into two three-dimensional subspaces, 
which are invariant under icosahedral symmetry opera- 
tions, the "physical" or "parallel" cutting space and 
the "perpendicular" space E 1 - . While the phonon displace- 
ment vector u is an element of E" , which transforms ac- 
cording to the icosahedral irreducible representation r 3 , 
the phason displacement vector w G E transforms ac- 
cording to the other three-dimensional icosahedral irre- 
ducible representation _T 3 . 

Both displacements only depend on the position vector 
x e EK The variations of this displacement fields are 
described by the strain fields 



u _ 1 / 9m duj 

11 2 I Ox, dx t 



dwj 



(1) 



In a linear theory the elastic energy is a quadratic form of 
the strain field: 

F(e) = -C a i/3je ai S0j a, (3=1... 6 i,j = 1...3. 

(2) 

Differentiation with respect to the strain e a i leads to Hooke's 
law 



dF 



(3) 
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where a denotes the generalized stress field. The meaning 
of the "phason stress" a w is not quite clear. According to 
the definition (||) it is the force conjugated to the phason 
strain e w . Phason strain is — in an atomic description — 
related to rearrangement of atoms on a small scale (so- 
called "flips"). Forces which inhibit such rearrangements 
might be considered as "pinning forces" . 

It is shown by group theoretical arguments M that 
Hooke's elastic tensor C for icosahedral symmetry con- 
tains only five independent elastic constants (/xi, ... , /is). 
From the transformation rules for u, x : _T 3 and w : r 
immediately follows jl3| the transformation behaviour of 
the strain fields 



i u : r 1 + r 5 



r 4 + r 5 



(4) 



which means that they can be written in an icosahedral 
basis 



„u5 



~w5 



(5) 



where e( u / w ) d is a d-dimensional vector transforming ac- 
cording to r d (see appendix A.l for details). In this sym- 
metry-adjusted coordinate system for the strains s and in 
an analogous way for the stresses a Hooke's law can be 
written as 



- a ul- 




a u5 








v w5 . 





"/ii 

/i 2 /i 3 

/i 4 

. /i 3 /i 5 



the Lame-constants /i, A by 

= 2/i + 3A [i 2 = 2(i, 





- £ ul- 




£ u5 








£ w5 . 






2 M , 



(0) 



(7) 



expressing the fact that an icosahedral quasicrystal be- 
haves like an isotropic medium when phason dynamics 
are frozen out. 

Hi and [i^ describe the elastic energy of pure phason 
strain. Contrary to the phonon case above, the phason 
elasticity is not isotropic, when these two constants are 
different. To simplify the calculations we assume in a first 
step that 



^4 = M5 



(8) 



is valid. We have made this assumption, because the SO(3) 
irreducible representation F l=i divides under icosahedral 
group operations into r 4 + T 5 , which corresponds to the 
splitting of the phason strain field e w (see eq. (f|)). If one 
demands not only icosahedral but spherical symmetry in 
the space of phason strains, the identity of the two pha- 
sonic elastic constants is necessary. Therefore we call eq. 
(§) "spherical approximation" . 

Phonon and phason elasticity are coupled by /i 3 . The 
assumption that this coupling is small, 



M3 < Ml>M2,M4 = M5, 



(9) 



allows us to develop a perturbation theory with respect to 

In ideal materials (without dislocations or defects) the 
only sources for stress fields are external volume forces f . 
Hence, the balance of forces is expressed by 



div cr + f = 0. 



(10) 



Please note, that the external force f — [f u , f w ], as the dis- 
placement [u, w], consists of two parts — the usual phonon 
force f" and the so-called phason force f w . Inserting Hooke's 
law (||) in the balance of stresses and forces (10) and ap- 
plying the definition (|l|) of the strain leads to generalized 
elastic equations: 



D(V) 



+ f = 0. 



(11) 



This system of six second order partial differential equa- 
tions describes the connection between applied forces and 
displacement field. In this paper we provide an approxi- 
mative solution by the method of Green's function. 

The differential operator D(V) is a 6 x 6- matrix and 
can be separated into four 3 x 3-matrices: 



D(V) 



D M ' U (V) D U ' W (V) 
D W < U (V) D w > w (\7) 



(12) 



For example, the second component of the phonon force is 
related to the displacements by D^iti + D^Wi + fa = 0. 
The phonon-phonon block D u,u is the well known differ- 
ential operator for the isotropic elastic continuum: 



D U,U (V) = [i A + {n + A) graddiv . 



(13) 



Applying the spherical approximation for the phason-pha- 
son-block explained above yields 



D W - W (V) = K\A, 
The phonon-phason coupling is described by 

D t0,Il (V) = (D"< t0 (V)) t 

F 1 (x,y,z) F 3 (x,y) F 2 (z,x) 



(14) 



A*3 

V6 



F 2 (x,y) Fi(y,z,x) F 3 (y,z) 
F 3 (z,x) F 2 (y,z) Fi(x,y,z) 



(15) 



with 



n , d 2 1 d 2 d 2 

F 1 (a,b,c) = ~— - -— + t— 



F 2 (a,b) 
F 3 (a,b) 



da 2 t db 2 
1 8 2 



dc 2 



2r 



t dadb 
2 



(16) 



dadb' 



2.2 Method of Green's function 

The Green's function method is commonly used to solve 
linear inhomogenous differential equations like (pi]). Ex- 
pressing the displacement field [u, w] (x) as a linear com- 
bination of the applied forces f , 



(x) = J G(x - x')f (x')dV, 



(17) 



J. Bachteler, H.-R. Trebin: Elastic Green's Function of Icosahedral Quasicrystals 



3 



leads to a system of differential equations for the "Green's 
function" G(x): 



D(V)G(x) + l«J(x) = 0. 



(18) 



The operator D and the Green's function G are 6x6- 
matrices, 1 is the six-dimensional identity operator and 
<5(x) the three-dimensional delta function. 

By Fourier methods this system of differential equa- 
tions can be transformed into a solvable system of al- 
gebraic equations p| . The problem is that until now a 
method for Fourier back transformation has not been pro- 
vided. 



2.3 Solution for spherical approximation and 
noninteraction 

If we assume that phonon elasticity is not coupled with 
phason elasticity (fa = 0), the differential equations ( [l8| ) 
become very simple: 



D u < u (V)Go'"(x) + 15(x) = 
D^(V)G"' m (x) + M(x) = 0. 



(19) 
(20) 



Now 1 is the three-dimensional identity operator and the 
Green's function matrix G is divided into four blocks by 
analogy with the differential operator D in (|l2|). The so- 
lutions of these two systems are given by 



GU.ll/ \ 
o ( x ) 



V) 



8 ,MAW (A + 3 ^H + (A+M) iy } 

(21) 
(22) 



where x ® x denotes the dyadic product of x with itself. 
Lord Kelvin [|| has provided solution (fj"i"l). Eq. ( gg ) is 
mathematically equivalent to the solution for the electric 
potential in a given charge density (see for example Feyn- 
man Q, II 6-1]). 

In the next section we develop a perturbation theory 
with respect to fa, where Go from (pil|), ( p2[ ) are the solu- 
tions of zeroth order. 



3 Solution for finite phonon-phason 
interactions 

Based on the solution for vanishing coupling constant (fa = 
0) given in the previous section we now develop an ap- 
proximative solution for small but fin ite c oupling between 
phonon and phason elasticity. Section 3J provides the per- 
turbation of the Green's function G(x) with r espe ct to fa 
and yields some relations for G(x). In section [T^ we con- 
struct an invariant basis for the perturbative solutions, 
which fulfills these relations. Finally, in section 3.3 we find 
the solutions up to the order 2 in fa. 



3.1 Perturbation theory 

The differential operator D(V) defined in ( |L?] ) consists of 
two parts 



D(V) = D (V) + Di(V) 



with 



D (V) 
Di(V) 



D«.«(V) 
D W > W (V) 

D U ^(V) 
D W ' U (V) 



(23) 

(24) 
(25) 



where the subscripts 0, 1 denote the order of fa (see (|15[)). 
Inserting (|23| ) and the expansion of the Green's function 
G 



G(x) = G (x) + Gi(x) + G 2 (x) 



(26) 



into the system of differential equations ( |18| ) and compar- 
ing the terms with the same order in fa leads to 



D (V)G (x) + H(x) =0 
D (V)G„(x) + Di(V)G n _i(x) = n 



1,2,. 



(27) 
(28) 



From the form (|2|) and of D and Di it follows im- 
mediately that the phonon-phason-coupling blocks of the 
Green's function of even order vanish, while the only re- 
maining ones are of odd order: 



Gu,w f~w ,11 f~u,u fW ,w ■ n 1 

2i - b 2i - b 2i+l - b 2i+l - U « - U, 1, . 



(29) 



Because equation (|27| ) of zeroth order is equivalent to the 
noninteracting equations ( |l9| ) (pp|), the nonvanishing part 
of the Green's function Go is given by (21) and (p2j). It is 
necessary to solve (28) to obtain the Green's functions of 
higher order. 

Let us now consider the analytic form of the solutions. 
Four properties result from the system of differential equa- 
tions (III): 



"Homogeneity" of degree —1 

G(ax) = a _1 G(x) aeM. 
"Inversion symmetry" 

G(-x) = G(x) 
"Transposing symmetry" 

Gjj(x) = Gji(x) 
"Icosahedral symmetry" 

G(x) = r 6 (g)G (r 3 (g)~ 1 x)r ,e (#) ~ 1 



(30) 



(31) 



(32) 



(33) 



The six-dimensional representation r 6 describes the trans- 
formation of the six-dimensional vector [u, w] . It is the 
sum of the two irreducible representations -T 3 and -T 3 . 
Each part G„ (x) of the perturbation series must have these 
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four properties, too. The relations ( pCj ) and ( |31| ) are valid 
for each matrix element G n -ij(x). Especially the matrix 
elements of the solution of zcroth order can be written as 



Go;y(x) 



1 



if E &r^y 

2=0 m=-l 



(34) 



where Y" ; m (x) are the "spherical harmonics" in cartesian 
coordinates. 

The solutions of higher order can be expressed in the 
same form with higher Z max . It is possible to make an 
ansatz with coefficients g^j and to transform (^8j) into 
an algebraic system of equations. But to avoid an explod- 
ing number of coefficients, we take into account that each 
solution G„ (x) must be invariant under icosahedral trans- 
formations (p3|). 

The separation of the Green's function G(x) into the 
four blocks G°'mx) with a, j3 S {u, w} transforms rela- 
tions (k32) and (B3l) into 



ij ~~ ," 



ij ~~ ji 



ij ~ ji 



(35) 



and 



G"'"(x) = r a (g)G^(r 3 (g)- 1 x)r^g)- 1 (36) 



with the definition F u := T 3 and r w := T 3 '. In the next 
section we construct icosahedral invariants for the four 
blocks, e.g. 3 x 3-matrix-functions, which fulfill relations 
©, ©, © and ©. 

3.2 Icosahedral invariants 

Any 3 x 3-matrix-function can be expressed as 

^M lSl (x) (37) 

i 

with 3 x 3-basis- matrices Mi and scalar functions s^(x). 
To achieve a basis in the space of icosahedral invariant 
matrix-functions it is necessary to divide the space of ma- 
trices and the space of scalar functions into icosahedral 
invariant subspaces, i.e. to construct the icosahedral irre- 
ducible bases (see for example |H| ) . By building the inner 
products of all possible combinations of irreducible bases 
of the matrices and the functions, where both bases belong 
to the same irreducible representation, we obtain a com- 
plete basis of the space of icosahedral invariant matrix- 
functions: 



J^(x) = 



1 



■E 

i=l 



(x) 



(38) 



In this formula a, (3 S {3,3'} specify the block and 7 
denotes the irreducible representation r 1 of dimension 
<i 7 , to which the irreducible bases {M"'' 3 ' 7 }^!. ^^ and 



{s'' 7 (x)}i = i...d T belong. The index n on the left hand side 



is a function of the representation index 7 and the or- 
bital quantum number I on the right. It is the enumer- 
ation of all possible (7, I) combinations, ordered first by 
I 6 {0,2,4,6,8} and second by 7 € {1,4,5}. In the next 
two paragraphs we discuss the two irreducible bases. 

The irreducible basis of the 3 x 3-matrices M can easily 
be calculated for all four blocks, but because of (|35|) it is 
sufficient to consider the three cases (u,w), (w,w) and 
(w,u). The "Clebsch-Gordan-Series" are given by 

- r 3 ® r 3 = r 1 © r 3 © r 5 

- r 3 ' <g> r 3 = r 4 © r 5 

- r 3 ' <g> r 3 ' = r 1 © r 3 ' © r 5 . 

The explicit form of the irreducible basis matrices is given 
in appendix A. 2. The matrices which belong to the ir- 



reducible representations r 3 and r 3 are antisymmetric. 
We do not need them for our purpose, because the Green's 
function must be symmetric (32. p5[). 



The space of homogeneous functions of degree con- 
sists of spherical invariant subspaces, denoted by the or- 
bital quantum number whose bases are given by Bi — 
{F i ' T1 (x)} TO= _;...;. The inversion symmetry ( pSl] ) allows us 
to limit the considerations to even quantum numbers I. 
The icosahedral group is a subgroup of the spherical group 
SO (3). The spherical irreducible basis B\ is icosahedrally 
reducible for each I > 2. We have calculated the icosa- 
hedral irreducible bases {s'' 7 (x)},...d„ from Bi for / = 



0, 2, 4, 6, 8. They are given in appendix A. 3 for I = 0, 2. 

Combining all possible pairs of irreducible matrices 
and irreducible functions using equation ( |38"1) leads to a set 
of icosahedral invariant matrix functions /"'^(x), which 
shows inversion ( pl| ) and transposing ( |35| ) symmetry. These 
invariants are homogeneous of degree 0. By dividing them 
by |x| one obtains homogeneity of degree —1 (|30|). 



3.3 Solution in terms of coefficients 

It is possible to express the Green's function of order in 
terms of the icosahedral invariant matrix functions calcu- 
lated in the previous section. Inserting 



l = ^/ 1 3 ' 3 (x) = % /3/f' 3 '(x) 
x®x 1 33 c.3,3, v 

in equations ( fH| ) and (S3) leads to 
1 



(39) 
(40) 



GU.lt/ \ 
W 



247r/i(A + 2/^)|x| 
{2^3(2A + 5m)/i' 3 (x) + 3V5(A + m)/ 2 3,3 (x)} (41) 



'(x) 



V3 



r 3',3' 



47rXi|x| 



(x). 



(42) 



The Green's functions of order 1 and 2 are calculated by 
inserting an ansatz in the invariant matrix functions into 
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the recursion equation ( pq ) and solving the systems of 
algebraic equations in the coefficients. The solutions are: 



Gil), 11/ \ 
i (x) 



2247rifi^(A + 2^)|x| 
{4V5(4A + f f/i)/i + 7\/l5(A + /i)/ 2 + 

10V21(X + n)h} (43) 



2 ( x ) 



Ml 



tt^ iM 2 (A + 2^) 2 |x| 



{"T^( 4a2 + 16 V + 19M 2 )A ~ tt(2A 2 + 8A^ + 5 M 2 )/ 2 



f26 
5\/21 
160f6 
5^ 



84 

(43A 2 + 242A^ + 342^ 2 )/ 3 + ^^(A + /i)(A + 3Ai)/4 



672 



1056 



(A + /i)(3A + 8/i)/ 5 



-(A + At ) 2 J 6 



262912 



35%/TT85 
7584 



(A + A1 ) 2 / 7 } (44) 



2 W ~ 7rJfM^ + 2/i)|x| 
{-y|(4A + H/i)/i - (A + 3„)J 2 + ^J(2A + 13 M )J 3 



25\/7 , T 25v / 55 / , 

H — (A + /x)l4H — (X + ii)I 5 \. (45) 

672 v w 1056 v w x v y 

/„ is an abbreviation for I"'^(pc) with appropriate a, (3 £ 
{u,w}. 



4 Results and discussion 

With the help of the Green's function calculated in the 
previous section it is now possible to examine the displace- 
ment fields u(x) and w(x) due to a force field. Inserting 
the point force 

f (x) = f <5(x) (46) 
into equation ( |l7j ) leads to the displacement field 



(x) = G(x)f . 



(47) 



Because of the homogeneity of degree —1 (30) the radial 
behaviour is given by 



( X ) OC 7^7 



(48) 



It is more interesting to examine the variation of the dis- 
placement fields on a sphere — the icosahedral symmetry 
should be reflected. Let us now consider a phonon force 



ftf. The corresponding displacement fields in perturbation 
theory of second order 



u(x) 
w(x) 



QU,U 
QW,U 



(x)f u 



Gr(x)f u + Gr(x)f u + ... (49) 
Gr tl (x)f " + ... (50) 



show, that the "isotropic" phonon displacement Gq'^x)!^ 
is corrected by a term of second order and that the phason 
displacement is given by a term of first order. We choose 
a spherical coordinate system with £ [0 . . . 2n] and 9 £ 
[0 . . . 7r] , which depends on the force direction {q in such 
a way, that fg always points to the "North Pole" (9 = 0). 

This force breaks the full icosahedral symmetry, but 
nevertheless the symmetry on parallels of latitude (9 = 
const) is preserved: n-fold symmetry, when the force points 
into such a direction of the icosahedron. The figures show 
the absolute value of the displacement fields u (figures [l] 
and ||) and w (figure ||) on such a parallel with 9 = j. 
For these figures we have choosen the elastic constants as 
given by Jaric and Mohanty Jl6| . 

All figures reflect the deviation from the spherical sym- 
metry but show the icosahedral symmetry: the two-fold 
axis in figure [l], five-fold in figure p] and three-fold in fig- 
ure |. In the case of phonon displacement this deviation is 
less than one percent. It is possible to examine the three 
components of the displacement fields or different values 
of 9, but the results are very similar. 

It is important to note that the deviation from the 
spherical symmetry is present although we have made the 
spherical approximation in the space of phason strains. 
Because the space of phonon strains is isotropic by defini- 
tion, this deviation is only due to the nonvashing coupling 
between these two spaces, depending in second order on 

We will now consider the five-fold axis (figure g) in 
detail. The absolute value of the phonon displacement field 
in second order on the parallel of latitude with 9 = \ is 
given by 



u 



(2) 



(0)|:=|Gr(r = l,0=7, 



)*o I 



MalQ 2 \l a °(-) + Qi(-)cos(50) + a 2 (-)cos(10(/>), 



A 4 



(51) 

where the coefficients ao, a±, and a 2 are functions in the 
ratio A//i. They are plotted in figure U as a function of 
Poisson's ratio v, which is connected to A and jj, by 



A 



2;/ 



l-2v 



(52) 



In the physically interesting interval v £ [0, 0.5] the coeffi- 
cient ao is half an order of magnitude larger than a\ , which 
again is one order of magnitude larger than a%. There- 
fore the maximum and the minimum of |u' 2 )(0)| are at 
the same positions as the maximum and the minimum of 
cos(50). The difference between maximum and minimum 
is then given by 



|u( 2 >(0)Mu( 2 )(^)|. 



(53) 
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In this paper we have discussed the theory of elasticity 
in icosahedral quasicrystals. We provide for the first time 
an approximative solution for the elastic Green's function, 
which reflects the deviation of the isotropy. For this ap- 
proximation we made the assumption, that 

— the two phason elastic constants are equal (^14 = /X5) 
and that 

— the coupling constant /13 is small. 

Measurements of Boudard, de Boissieu et al 
in AlPdMn show that the first assumption is not fulfilled. 
The constant K 2 should be zero (see appendix |b|). Never- 
theless, the qualitative results of this paper remain valid, 
because they are based on the coupling between phonons 
and phasons. The component of the phason strain e w5 , 
which is coupled to the phonon strain, is directly con- 
nected to the constant /X5. The value of the other phason 
constant does only indirectly affect this coupling. 

The elastic Green's function is provided up to second 
order of perturbation theory in a basis of icosahedral in- 
variants. With the help of this Green's function we exam- 
ined the phonon and phason displacement fields due to a 
phonon point force. Although the icosahedral symmetry 
is close to the isotropic symmetry, the deviation from the 
isotropy is observable. 

A Icosahedral basis 

A.l Icosahedral irreducible strain 



We have choosen the same coordinate system as in 
[l9|| . Using the golden number r = (1 + v / 5)/2 the 
hedral irreducible strain is given by: 

' 1 f(£n+£22 + £33) 



e = 



1 

7s 



1 

if 21 

32 



-33 



-c 

Is 



•reS 



■ re 



~w5 



1 

76 



f (-r 2 ^ 4 


1 -U 

T 2 t 22 ' 


2 W E H 


~re u 22 




V2£i2 














F w _|_ p w 
fc ll ' fc 22 


Zfc 33 


>/2(re3i - 


T £ 12> 


V2(ts$ 2 - 


~ 


_V2(rer 3 - 





icosa- 
(54) 

(55) 
(56) 



(57) 



A. 2 Icosahedral irreducible matrices 

In order to write the icosahedral irreducible matrices in a 
compact form it is useful to define a "permutation opera- 
tor" 



P:Ih P(M) with P{M) i;j := Mi_i,j_i 



mod 3 



. (58) 



which permutes the matrix indices modulo 3. 
— phonon-phonon-block: 



M- 



3,3;1 



V3 
3 



3.3:3 



AL 



M 3,3;5 



10 

1 
1 

3,3;3\ 



AL 



3,3:3 



V2 
2 



AL 



3,3:3 




1 

0-10 

3,3;3\ 



p2 ( - M 3,3; 



AL. 



AL 



3,3;5 
3,3:5 



V3 
6 

1 

2 



— T U 

£ 

Or 

i 0" 
-r 
1 

3.3:5\ 



P(M 3 , 
phason-phonon-block 



M 3,3;5 = V2 

3 2 

3,3:5 



10 

1 




AL, 



P 2 (M 3 ' 3;5 ) 



71/ 



3',3;4 



AL 



3',3;4 



M 3,3;l 



P(M 2 3 '' 3;4 ) 



M 3 '< 3;4 = ^ 



V3 
3 



t 




AL, 



3'. 3:4 



p2 (M 3',3;4 ) 



M 



3',3;5 



AL 



3',3;5 



V2 
2 

V3 
3 



AL: 



3',3;5 



-10 

1 
00 

i 0' 

-t0 


= P 2 (Af 3 3 '' 3;5 ) 



AL 



3', 3:5 



AL, 



3',3;5 



V6 
6 

P(AL 



-10 
0-10 
2 



3 ,3; 5) 



The icosahedral irreducible matrices of the phason- 
phason-block can be generated from the matrices of 
the phonon-phonon-block by the substitution 



r 4 



p d 



(59) 



All matrices are normalized with respect to the quadratic 
vector norm, where {<%|i, j = 1 ... 3} is considered to be 
the normalized basis of the vector space. 



A. 3 Icosahedral irreducible functions 



s{ (x) = 1 



V3 



S 2;5 (x) = -^[-r 2 .* 2 + + (r + -)z 2 ] 

DX Z T z T 

2;5/ \ 1 r -L 2 2 1 2i 

s 2 ( x ) = -^\- x ~ T y +z \ 



(60) 

(61) 
(62) 



2:5/ \ V2 2 ;5/ V% 2;5/ V% 

S 3 ( X ) = — ^ S 4 ( x ) = — V Z S 5 ( X ) = — zx 



(63) 
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All functions are normalized with respect to the quadratic 
vector norm, where {\J - v x i z k > 0,i + j + k = 
1} is considered to be the normalized basis of the vector 
space, which is the space of all spherical harmonics to the 
orbital quantum number I. 



B Elastic constants 

Unfortunately, the definitions of the elastic constants of 
icosahedral quasicrystals differ in the literature. Therefore 
we want to compare our definition given in ([|) to others. 

- Trebin, Fink, Stark @|3 



Mi = 3Ai 

M4 = ^4 



M2 = A 2 

M5 = ^5 



M3 



(64) 



- Shaw, Elser, Henley pl| , Widom |gj and Henley 

fuse the well known Lame-constants A and n, see 
The three remaining constants are: 



M3 



= K 3 V6 



Ms = Ki ~ 

(65) 



Boudard, de Boissieu et al JT?] . [T^l use the same 
elastic constants (see reference 19 in ||18|]). 
Ding et al ylj also use the Lame-constants, but the 
"phasonic" constants are different: 



M3 



Ry/6 ( i i =K 1 - 2K 2 Ms = K 1 



K 2 . 

(66) 
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Fig. 1. Absolute value of phonon displacement |u| versus <f> for 
9= j, 2-fold axis. 
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Fig. 2. Same as figure |l[ but 5-fold axis. 
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Fig. 3. Absolute value of phason displacement |w| versus 4> for 
9 = j , 3- fold axis. 
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Fig. 4. Coefficients a ,ai, and as a function of Poisson's 
ratio v. 



